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Effects  of  Spatial  Light  Modulator  Transmissive 
and  Reflective  Dead  Zones  on  Optical  Correlation 


1.  INTRODUCTION 


We  previously  analyzed  the  effects  of  opaque  dead  zones  (ODZs)  in  pixels  of  spatial  light 
modulators  (SLMs)  used  in  standard  4-f  or  2-f  optical  correlators.  ^  There,  we  showed  by 
computer  simulation  and  approximate  analytical  calculations  how  changes  in  these  ODZs 
affected  the  correlator's  signal-to-noise  ratio,  correlation  peak  intensity,  and  the  efficiency 
of  the  system's  light  throughput  when  either  a  phase-only  filter  (POP),  a  binary  phase-only 
filter  (BPOF),  or  a  classical  matched  filter  (CMP)  is  used  in  the  filter  plane. 

Our  goal  here  is  to  examine  how  these  same  correlator  performance  characteristics  are 
affected  when  pixels  in  the  dead  zone  portions  of  the  SLM,  rather  than  completely  blocking 
the  incident  light,  completely  transmit  (or  reflect)  it  towards  the  output  detector.  Since  the 
mathematics  describing  the  reflecting  mode  is  identical  to  that  describing  the  transmitting 
mode,  we  will  refer  to  both  of  them  as  transmissive  dead  zone  (TDZ).  In  this  report  we 
provide  accurate  analytical  calculations  and  additional  simulations  for  TDZ  using  the  same 
three  filters  as  above.  We  will  also  show  that  the  analytical  treatment  of  TDZ  is  a  general 
one  because  it  includes  the  ODZ,  as  well  as  the  mixture  of  ODZ  and  TDZ,  as  special  cases. 

Received  for  Publication  13  April  1992 

1  Gianino,  P.D.  and  Woods,  C.L.  (1992)  Effects  of  spatial  light  modulator  opaque  dead 
zones  on  optical  correlation,  Appl  Opt,  31:4025-4033. 
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Furthermore,  we  allow  for  the  possibility  of  placing  a  dc  block  at  the  center  of  the  filter 
plane  SLM.  For  the  simulation  we  follow  the  same  sequence  of  computer  operations  shown 
In  the  computer  flow  diagram  in  Figure  2  of  Reference  1 . 


2.  BACKGROUND 


The  SLMs  in  the  input  and  filter  planes,  designated  as  SLMl  and  SLM2,  respectively,  and 
the  detector  in  the  output  plane  are  assumed  to  have  arrays  of  64x64  pixels.  We  adopt  the 
model  that  each  pixel  of  a  SLM  is  subdivided  into  two  regions.  One  region,  called  the  "active 
zone"  (AZ),  encodes  the  addressed  amplitude  or  phase  signal  modulating  any  light  incident 
on  it.  The  remaining  region  of  the  pixel,  called  the  "dead  zone"  (DZ),  does  not  encode  any  of 
the  addressed  cimplltude  or  phase  signal.  In  our  model  this  dead  zone  is  either  opaque 
(blocking  100  percent  of  any  light  incident  upon  it),  or  transparent  (transmitting  or 
reflecting  all  of  this  light).  Examples  of  SLMs  whose  dead  zones  are  opaque  are  SEED, 
symmetric  SEED,  magneto-optic,  Fabry-Perot  etalon  and  PLZT  devices.  Liquid  crystal  TVs 
have  transmissive  DZs,  while  the  deformable  mirror  devices  have  reflective  DZs.  Other 
devices,  such  as  liquid  crystal  light  valves,  microchannel  SLMs  and  ferroelectric  liquid 
crystal  SLMs,  have  more  complicated  DZs.  but  these  could  be  approximately  modeled  by 
assuming  an  average  value  of  the  complex  amplitude  transmittance. 


3.  MATHEMATICAL  ANALYSIS 


3.1  Input  Plane 


A  diagram  of  a  typical  pixel  in  SLMl  is  shown  in  Figure  1:  Xj,  yj,  are  the  spatial 
coordinates  in  the  input  plane.  The  AZ  (the  light  area)  is  positioned  at  the  center,  while  the 
transmissive  (or  reflective)  DZ  (the  dark  area)  is  located  around  the  edge.  The  dimensions  of 
the  pixel  are  given  by  Cj^,  Cy;  those  of  the  AZ  by  b^^.  by.  A  similar  diagram  describes  the 
typical  pixel  in  the  filter  plane,  except  that  the  dimensions  of  the  pixel  are  given  by  h^,  hy 
and  those  of  its  AZ  by  e^.  Cy.  The  fractional  areas  for  the  dead  zones  in  both  SLMs  are 
specified  as 


Zj  =  1  -  b^by/Cj^Cy  (for  SLMl)  lU 

Z2  =  1  -  e^ey/h^hy  (for  SLM2)  (2) 
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Cx 


Figure  1 .  Diagram  of  a  Typical  Pixel  of  Dimensions  Cj^,  Cy  in  SLM 1 .  The  active 
zone  in  the  center  has  dimensions  b^.  by.  The  darkened  area  corresponds  to  the 
transmissive  (or  reflective)  dead  zone. 


The  spatial  coordinates  in  the  filter  plane  are  given  by  X2.  y2.  while  the  spatial  frequency 
coordinates  in  that  same  plane  are  given  by  ^2-  'Hz- 

A  block  diagram  depicting  the  various  system  components  and  signals  involved  in  the 
correlation  process  is  shown  in  Figure  2;  sCx^,  yj)  is  a  continuous  input  signal.  The  final 
signal  amplitude  encoded  on  SLMl  due  to  the  active  areas  in  the  pixels  is  the  function 
w(xj,  yj^).  According  to  our  model,  this  function  is  constant  within  a  given  active  zone  (the 
constant  value  representing  a  sampled  version  of  s)  and  is  zero  in  the  DZ  (a  detailed 
mathematical  description  of  w  and  the  other  functions  to  be  introduced  will  be  given  in  the 
Appendix).  The  DZ  of  each  pixel  in  SLMl  transmits  (or  reflects)  100  percent  of  the  light 
illuminating  it.  This  light  from  the  DZ  transmittance  is  represented  by  tglxj.  yj).  The  sum 
of  the  exiting  signals  w  and  tg  are  Fourier  transformed  (FT)  by  lens  1  into  the  sum  of 
transforms  W(^2>  Ts(^2-  'Hz^  respectively,  where  (^2-  Tz)  spatial  frequency 

coordinates  in  the  filter  plane. 
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Figure  2.  Block  Diagram  of  Correlator  Showing  the  Different  Components 
and  Signals  Involved  in  the  Correlation  Process.  The  various  signals  are 
explained  in  the  text. 


3.2  Filter  Plane 

The  function  F*(^2-  '•12)  Is  the  continuous  filter  transform  function  for  one  of  the  three 
desired  filters.  The  final  amplitude  of  the  filter  transform  function  encoded  on  SLM2  due  to 
the  active  areas  in  the  pixels  is  the  function  G(^2*  ^2^-  According  to  our  model,  this  function 
is  constant  within  a  given  AZ  (the  constant  value  representing  a  sampled  version  of  F*)  and 
is  zero  in  the  DZ.  Analytically.  F*  is  formed  from  the  Fourier  transform  of  s(xi.  yj)  and  G  is 
obtained  from  F*  (see  Eq.  (A- 12)  in  the  Appendix). 

The  DZ  of  each  pixel  in  SLM2  not  only  reduces  the  optical  efficiency  of  the  filter,  but  it 
also  Introduces  noise  by  transmitting  (or  reflecting)  all  of  the  incident  light.  This  unwanted 
DZ  transmittance  is  designated  by  Tp(^2>  ^2^  Furthermore,  a  transmissive  dead  zone  in  a 
SLM  acts  like  a  grating,  transmitting  a  significant  part  of  its  energy  into  the  (0.  0)  order  and 
diffracting  the  rest  into  higher  order  images.  The  physical  size  of  SLM2  is  chosen  so  as  to 
transmit  only  the  (0,  0)  orders  of  W  and  T^  onto  SLM2.  These  transmitted  portions  are 
designated  by  W  and  Ts,  respectively.  The  remaining  diffraction  orders  fall  outside  of  the 
SLM  and  are  blocked.  Since  the  Tg  portion  contributes  noise  and  is  localized  around  the 
optic  axis,  an  opaque  dc  block  (at  the  (0,  0)  frequency  position)  eliminates  a  great  deal  of 
this  input  TDZ  noise  from  the  system.  This  dc  block,  whose  transmittance  is  given  by 
Tdc(^2>  'Hz)*  is  indicated  by  the  small  dark  rectangle  at  the  center  of  SLM2  in  Figure  2.  The 
net  signal  exiting  from  SLM2  consists  of  four  product  terms:  a  signal  term  W'.G.Tdp.  two 
imaged  noise  terms  (W'+Ts).Tp.Td(.,  and  a  filtered  noise  term  Ts.G.T^c.  The  resulting  four 
terms  are  transformed  by  lens  2  to  produce  a  total  correlation  amplitude  pattern  kj(x3,  y3) 
at  the  detector,  where  X3,  y3  are  spatial  coordinates  in  the  detector  (or  output)  plane. 
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3.3  Output  Plane 


The  total  correlation  amplitude  k-p  can  be  written  as  the  sum  of  the  inverse  Fourier 
transforms  of  the  four  product  terms  mentioned  above.  These  four  terms  may  be  written 
explicitly  as 


kc(x3.  yg)  H  FTMW'.G.TdJ  =  [w'(x3.  y3)  *  g(x3.  y3)l  O  tdc{x3.  y3) 

k3(X3.  y3)  =  FTMW'.Tp.Tj^l  =  [w'{x3.  y3)  O  tp(X3.  y3)l  O  t,,e(x3.  y3) 
kR(x3.  y3)  =  FT-MTs.Tp.Tdcl  =  lts(x3.  Ya)  ©  tp{x3.  y3)l  O  t,,^(x3.  yg) 
kpRtxa.  Ya)  =  FrMTs.G.T,,J  =  [ts(x3.  y3)  •  g(x3.  y^)]  O  t„^(x3.  y3) 


where  FT"!  refers  to  the  Inverse  Fourier  transform  operation.  •  represents  the  correlation 
operation,  O  the  convolution  operation  and  w'.  ts.  g,  tp  and  t^jp  are  the  inverse  Fourier 
transforms  of  W.  Tg.  G.  Tp  and  T,,p.  respectivelY-  Equations  (3)  and  (6)  contain  the 
correlation  operation  because  of  the  complex  conjugate  term  in  the  g  function. 

For  these  and  subsequent  equations  we  employ  the  following  notation  for  certain 
functions  that  occur  repeatedly  in  both  the  input  and  detector  planes: 


v„  =  rect(x/bv)/bx 


=  sinc(x/Cx)/Cx 

mjj  =  sinc(exX/hxLx) 

r^  s  rect(x/Lx) 


(7) 


(8) 


(9) 


(10) 


in  which 
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sinc(a)  s  (sinmj/jta. 


(11a) 


and 


reel 


0.  |x|  >  Xq 

1.M<Xq 


(11b) 


We  also  utilize  the  fact  that  Qi  O  Qi  =  Qi  for  l=x  or  y.  Is  the  x  dimension  length  of  SLMl  as 
well  as  of  the  detector,  and  the  variable  x  can  refer  to  either  xj  or  Xg.  depending  upon  which 
plane  Is  being  referred  to.  There  will  be  another  set.  namely.  Vy.  Qy.  my.  and  ry,  which  have 
the  same  functional  forms  as  above,  with  all  x  subscripts  and  variables  replaced  by  y.  In 
the  limit  of  small  pixel  and  active  zone  sizes  (that  is.  b.  c.  e.  h«L  and  e«h).  functions  v  and 
q  of  Eqs.  (7)  and  (8)  become  8  functions,  and  functions  m  and  r  of  Eqs.  (9)  and  (10)  become 
flat  over  the  SLM  faces. 

The  explicit  2D  forms  for  the  functions  w'.  tg.  g.  tp,  and  t^^.  which  are  mentioned  in  Eqs, 
(3)-(6).  as  well  as  their  Fourier  transforms,  are  derived  in  the  Appendix.  Their  final  forms 
are  given  in  Eqs.  (A-7).  (11).  (13).  (17).  and  (20).  respectively.  When  these  are  substituted  into 
Eqs.  (3)-(6)  and  only  the  zero  diffraction  order  terms  are  kept,  there  results 


kc(x3.  Ya)  =  (1-Zi)  (I-Z2) 


(12) 


{[(^%  yJ  ®  Vy)  •  Vy]  ■  [■■’ (-V  -^3)  ■  ”’x“y]  ®  lx",  }  ®  ‘do 


,)  ^2  T[(*(*3-  ^3)  ®  Vy)  ■  Vy]  ®  <!x'ly  J-  ®  ‘do 


■Zj-'*'' (*3-^3)® ‘do 


ya)  =  Z.Z2  [  Vy  ®  'Ix'Iy]  ®  ■  Z2  ■  )  (%  '‘3)  ® 

ya)  =  Z,(l-Z2)  {  (v^  '[f"  (-X3.  -y,)  •  ni.mj)  O  q,q,}  O 


(13) 


(14) 


(15) 
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The  Inverse  Fourier  transform  of  F*.  denoted  by  f*.  is  also  known  as  the  impulse  response  of 
F*  or  the  filter  reference  image. 

These  equations  reveal  the  basic  operation  of  the  system.  They  clearly  show  that  the 
final  correlation  consists  of  a  "true"  correlation  term  (the  one  containing  s»n,  plus  three 
background  terms.  The  true  correlation  term,  Iq,  of  Eq.  (12),  is  produced  solely  from  the 
active  zones  in  the  input  and  filter  SLMs.  It  is  identical  to  the  ODZ  correlation  discussed  in 
Reference  1  and  given  specifically  by  Eq.  (6)  there.  Under  these  conditions,  Eqs.  (4)-(6)  and, 
therefore,  Eqs.  (13)-(15),  vanish  because  t3(xi,  yi)=Tp(^2-  The  three  remaining  k  terms 

in  Eqs.  (13)-(15)  contribute  to  the  background  noise  and  coherently  interfere  with  the  k^ 
signal  term.  For  example,  the  "signal"  image  background  term  kg  is  the  image  of  the  input 
signal  that  has  passed  through  the  DZ  area  of  the  filter  SLM.  The  appearance  of  input 
images  in  the  output  plane  has  been  noted  by  Goodman^  for  the  VcUider  Lugt  correlator  and 
by  Flavin  and  Horner^  for  the  amplitude-encoded  phase-only  filter.  Neither  of  these  cases, 
however,  exhibit  the  specific  TDZ  behavior  considered  here.  The  "rectangular  input  TDZ 
background  term  kj^  is  the  image  of  the  rectangular-shaped  TDZof  the  input  SLM  that  has 
also  passed  through  the  TDZ  area  of  the  filter  SLM.  The  "filtered  rectangular  input  TDZ  " 
background  term  kpp  is  the  light  from  the  TDZ  of  all  pixels  in  the  input  SLM  that  has 
passed  through  the  filter's  AZ  areas.  The  energy  of  this  light  is  concentrated  near  the  on- 
axls  pixels  of  the  filter,  producing  a  slowly  varying  pattern  in  the  correlation  plane.  As 
indicated  in  the  above  equations,  these  three  noise  terms  are  proportional  to  a  product  of  a 
TDZ  and  either  an  AZ  or  a  TDZ. 

Because  of  the  pixellation  of  the  SLMs.  all  4  terms  have  replicated  images,  with  their 
(0,  0)  order  on  axis  and  weaker  orders  off  axis.  For  example,  kp  is  a  series  of  2D  rect 
functions  in  the  output  plane,  with  the  central  one  being  the  highest  and  the  others 
diminishing  in  height  as  they  get  farther  away  from  the  center.  The  kg,  kp,  and  kpp  terms 
coherently  interfere  with  and  increase  the  correlation  peak  and  may  also  contribute 
secondary  maxima.  Such  higher  order  peaks  were  detected  in  computer  simulations  by  Yu'^ 
and  in  experiments  by  Javidi.®  In  addition,  both  have  been  detected  in  experiments  and 
analyzed  mathematically  by  Davis,  et  al.®  All  of  these  investigations  involved  systems  that 
employed  a  binarized  version  of  the  joint  transform  correlator  (JTC).  With  a  dc  block,  kg, 
kp,  and  kpp  are  high-pass  filtered  (or  edge-enhanced),  thus  decreasing  the  noise. 


2  Goodman,  J.W.  (1968)  Introduction  to  Fourier  Optics.  McGraw-Hill  Book  Co., 

NY,  Chap.  7,  pp.  174-177. 

3  Flavin.  M.A.  and  Horner.  J.L.  (1989)  Amplitude  encoded  phase-only  filters.  Appl. 

Opt,  28:1692-1696. 

4  Yu,  F.T.S..  Cheng.  F..  Nagata,  T.  and  Gregory.  D.A.  (1989)  Effects  of  fringe  binarization 
of  multiobject  joint  transform  correlation,  Appl.  Opt,  28:2988-2990. 

5  Javidi.  B.  and  Horner.  J.L.  (1989)  Multifunction  nonlinear  signal  processor: 
deconvolution  and  correlation.  Opt  Eng.,  28(8):837-843. 

6  Davis.  J.A..  Merrill,  E.A..  Cottrell,  D.M.  and  Bunch.  R.M.  (1990)  Effects  of  sampling  and 
binarization  in  the  output  of  the  joint  Fourier  transform  correlator.  Opt  Eng.,  29(9):  1094- 
1100. 
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In  the  Appendix  we  show  that  the  functions  appearing  in  the  filter  plane  are  limited  by 
the  finite  size  of  SLM2  and  are,  therefore,  low-pass  filtered  in  this  plane.  This  low-pass 
filtering  produces  the  q  terms  in  Eqs.  (12)-(15).  Note  that  the  filter  reference  image  position 
changes  the  correlation  peak,  as  described  in  Eq.  (12),  by  the  factors  m^-my.  The  qualitative 
results  that  can  be  gleaned  from  Eqs.  (12)-(15)  are  analyzed  more  fully  in  Section  5. 


4.  COMPUTER  SIMULATIONS 


For  computational  purposes  each  SLM  pixel  is  simulated  by  sampling  at  16  points, 
arranged  in  a  square  array  having  4  points  on  a  side.  Of  these  16  points  a  certain  fraction 
Zj  is  considered  to  make  up  the  dead  zone  (DZ)  area  of  SLMj,  and  the  remaining  fraction 
(1-Zi)  is  considered  to  make  up  the  (AZ)  area  of  SLMj.  The  4x4  point  sampling  of  the  64x64 
pixels  in  both  input  and  filter  planes  produces  arrays  of  256x256  points.  To  simultaneously 
produce  this  resolution  in  both  planes,  additional  zeroes  must  be  added  to  bring  the  final 
eurays  to  1024x1024  points.^ 

The  computer  simulations  used  the  same  input  signal  as  in  Reference  1  (namely,  a  girl's 
face)  with  each  of  the  three  filters.  The  quantities  computed  from  the  resulting  correlation 
pattern  were  its  peak  intensity  value  and  its  signal-to-noise  ratio  (SNR).  The  percentages  of 
the  energy  incident  on  or  lost  at  the  three  different  planes  of  the  system  were  also  evaluated. 
Three  different  cases  involving  TDZ  were  studied  and  are  listed  in  Table  1 :  Case  A,  in  which 
both  SLMs  have  equal  amounts  of  DZ;  Case  B.  in  which  only  SLM  1  has  DZ;  and  Case  C,  in 
which  only  SLM2  has  DZ.  These  cases  are  the  same  as  in  Reference  1  except  that  in  this 
study  the  number  of  sampling  points  per  pixel  eire  always  fixed  at  4  on  a  side  (Ni=N2=4). 


Table  1.  TDZ  Cases  Studied 


CASE 

SLMl 

SLM2 

A 

Zj  varies 

II 

B 

Zj  varies 

Z2  =  0 

HnsiiMl 

C 

Z,  =0 

Z2  varies 

If  Z,  =  Z2  =  0  for  emy  case 

K 
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It  should  be  mentioned  that  other  investigators^  ®-^  have  set  up  mathematical  models  to 
compute  how  pixellation  and  dead  zones  in  real  SLMs  such  as  deformable  mirror  devices 
and  magneto-optic  SLMs  affect  correlation.  Their  systems  were  of  both  the  4f  and  Joint 
Transform  Correlator  types.  In  each  case,  however,  their  computer  simulations’’  ®  ®  and 
experiments®  were  limited  to  that  one  particular  configuration  of  DZ  peculiar  to  each  type  of 
SLM  and  presented  no  careful  analytical  correlation  calculations.  No  investigations  have 
studied  the  effects  due  to  changes  in  DZ. 


4.1  Correlation  Peak  Intensity 


Figure  3  shows  how  the  computer-calculated  normalized  correlation  peak  intensity 
lp(=k^  changes  with  equal  input  and  filter  TDZ  for  Case  A  using  each  of  the  three  filters. 

The  abscissa  Zj  refers  to  the  TDZ  in  both  SLMs  since  Zj=Z2=Zj  for  Case  A.  These  intensities 
are  normalized  to  that  value,  neimely  (Ip=8.31)  for  the  POF  when  Zj=Z2=0  and  no  dc  block  is 
Inserted  in  the  filter  plane. 

The  solid  lines  in  Figure  3  give  the  Ip  for  the  ODZ  condition  without  a  dc  block.  These 
are  the  same  as  the  three  solid  curves  in  Figure  6  of  Reference  1  and  are  repeated  here  for 
reference.  A  dc  block  is  not  required  for  ODZ  since  no  noise  is  present.  The  dashed  lines 
give  the  Ip  for  the  TDZ  condition  with  a  1 -pixel  dc  block  placed  at  the  center  of  SLM2.  The 
dashed  lines  are  generally  higher  than  the  solid  lines  because  the  additional  light  getting 
through  the  transparent  dead  zones  Increases  Ip.  For  Zj=0  the  dc  block  significantly 
decreases  only  the  CMF  correlation  peak.  For  both  of  the  phase  filters  the  dotted  and  solid 
lines  are  fairly  close  up  to  DZs  of  about  20  percent.  Just  as  in  the  opaque  condition,  the  Ip 
decrease  as  TDZ  increases.  However,  at  some  point  the  Ip  go  through  a  minimum  and  begin 
to  increase.  This  occurs  at  DZs  of  about  20  percent  for  the  CMF  and  about  75  percent  for 
both  of  the  phase  filters.  This  happens  because  large  amounts  of  light  (which  is  really 
noise)  from  the  TDZs  are  added  to  the  correlation  pattern.  Finally,  when  TDZ  reaches  100 
percent,  all  three  dashed  curves  come  together  at  a  normalized  intensity  of  0.07.  This  seems 
reasonable  because  at  100  percent  TDZ  both  SLMs  are  completely  transparent:  transmitting 
all  of  the  incident  light  (represented  by  tg)  but  no  signal,  and  rendering  any  filter  function 
meaningless. 


7  Cravatt,  J.E.  and  Giles.  M.K.  (1990)  An  improved  model  of  the  TI  deformable  mirror 
device,  Proc.  S.P.I.E.,  1296:110-124.  Advances  in  Optical  Information  Processing  IV,  D.R.  Pape 
(ed.). 

®  Giles.  M..  Taylor.  J..  Grijalva.  N.  and  Gioannini.  B.  (1987)  Optical  image  correlation 
using  a  deformable  mirror  device:  a  feasibility  study.  Proc.  S.P.I.E.,  753:72-81,  Acousto- 
Optic.  Electro-Optic  and  Magneto-Optic  Devices  and  Applications,  S.  Lucero  (ed.). 

®  Bock,  B.D.,  Crow.  T.A.  and  Giles.  M.K.  (1990)  Design  considerations  for  miniature  optical 
correlation  systems  that  use  pixelated  input  and  filter  transducers,  Proc.  S.P.I.E.,  1347:297- 
309,  Optical  Information-Processing  Systems  and  Architectures  II.  B.  Javidi  (ed.). 
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Figure  3.  Normalized  Correlation  Peak  Intensity  !„  (Dashed  Lines) 
Versus  Transmissive  Dead  Zone  Area  (in  Percent)  lor  Case  A  With  a 
1-Pixel  dc  Block  Inserted  in  the  Filter  Plane.  The  calculated  values 
are  indicated  by  black  squares.  The  solid  curves  are  for  the  corre¬ 
sponding  opaque  dead  zone  mode  with  no  dc  block  included.  The 
dotted  curves  give  the  normalized  maximum  intensities  of  those  sub¬ 
sidiary  peaks  higher  than  Ip.  Their  calculated  values  are  indicated 
by  black  circles. 
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The  noise  from  the  dead  zones  also  enhances  secondary  peaks  in  the  correlation  pattern 
and  for  sufficiently  large  TDZs  can  even  cause  these  secondary  intensities  to  exceed  the  Ip. 
The  intensity  of  the  largest  of  these  secondary  peaks,  given  by  l,„jjx'  represented  by  a 
dotted  curve  in  Figure  3  for  each  filter.  Here,  the  dotted  l,„jjx  curve  for  the  CMF  begins  to  rise 
above  the  Ip  value  at  a  TDZ  =  20  percent.  The  corresponding  values  for  the  BPOF  and  POF 
are  approximately  65  percent  and  75  percent,  respectively.  Beyond  these  values,  the  dotted 
curves  rise  very  rapidly.  At  100  percent  TDZ  we  see  from  Eqs.  (12)-(15)  that  kj  reduces  to  kj^ 
only.  Hence,  the  secondary  peak  is  simply  the  meiximum  of  the  dc  block-filtered  image  of 
the  transparent  input  SLM,  which  in  the  present  situation  is  0.175. 

Similar  information  for  TDZ  in  the  input  plane  only  (Case  B)  is  shown  in  Figure  4.  The 
solid  curves  come  from  the  Case  B  ODZ  curves.  ^  The  dashed  curves  representing  the  TDZ  Ip 
do  not  experience  a  minimum  but  monotonically  decrease  in  magnitude  with  Increasing 
TDZ.  10  For  both  of  the  phase  filters  the  dashed  curves  stay  fairly  close  to  the  solid  curves  up 
to  D2^  of  about  50  percent.  Since  Z2=0  in  Case  B,  at  100  percent  TDZ  (Zj=l),  Eqs.  (12)-(15) 
show  that  kj  reduces  to  kp^  only.  Using  the  form  of  kpR  shown  in  Eq.  (6),  we  see  that  kp^  is 
proportional  to  (constant  •  g)  =  (constant). G(O.O).  Using  Eq.  (A- 12)  in  the  appendix  this  term 
eventually  reduces  to  (constant). F’(0,  0),  where  G(0,  0)  and  F*(0,  0)  are  evaluated  at  ^2”h2“0- 
Consequently,  at  100  percent  TDZ  each  filter  produces  its  own  characteristic  Ip  and  I^jax- 
Since  the  CMF  pixels  are  attenuating,  its  Ip  (and  l.nax)  1®  the  lowest  of  the  three  filters.  On 
the  other  hand,  the  POF  pixels  do  not  attenuate  so  its  Ip  (and  l,„ax)  is  the  highest  of  the  three. 
Although  the  BPOF  pixels  don't  attenuate  either,  they  divide  the  energy  between  a 
convolution  and  a  correlation  pattern  so  that  its  value  of  Ip  (and  l^ax^  ^®  lower  than  that  for 

the  POF. 

Similar  Information  for  TDZ  in  the  filter  plane  only  (Case  C)  is  shown  in  Figure  5.  The 
solid  curves  come  from  the  Case  C  ODZ  curves.  >  Just  as  in  Case  B  the  dashed  curves  decrease 
monotonically,  only  here  they  are  much  closer  to  the  solid  curves.  This  indicates  that 
having  TDZ  only  in  the  filter  plane  introduces  less  noise  into  the  system  than  having  TDZ 
only  in  the  input  plane,  as  in  Cases  A  and  B.  In  fact,  for  both  of  the  phase  filters  the  two 
curves  stay  fairly  close  together  over  the  entire  range.  Since  Zi=0  in  Case  C,  at  100  percent 
TDZ  (Z2=l)  then  kT-=ks  only,  which  from  the  form  shown  in  Eq.  (4)  is  proportional  to  w'. 
Thus,  the  Ip  and  I^j^x  are  equal  to  the  corresponding  intensities  on  a  dc  block-filtered 
version  of  the  signal  image. 


Gianino,  P.D.,  Homer,  J.L.  and  Woods.  C.L.  (1990)  Effects  of  SLM  transmissive  dead 
zones  on  optical  correlation,  Proc.  S.P.I.E.,  1347:240-246.  Optical  Information-Processing 
Systems  and  Architectures  II,  B.  Javidi  (ed.). 
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TRANSMISSIVE  DEAD  ZONE.  Z2  (%) 


Figure  4.  Same  as  Figure  3  Except  for  Case  B. 
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1.0 


TRANSMISSIVE  DEAD  ZONE.  (%) 


Figure  5.  Same  as  Figure  3  Except  for  Case  C. 
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As  a  performance  guide  for  the  two  phase  filters,  we  indicate  those  areas  of  TDZ  that  will 
cause  Ip  to  decrease  by  6  db  (that  is.  by  a  factor  of  4)  below  its  zero  dead  zone  value.  These 
operating  points  are  labelled  by  the  large  black  dots  in  Figures  3,  4  and  5.  This  condition 
occurs  at  a  TDZ  value  of  approximately  39  percent  for  Case  A  (in  which  b^  SLMs  have 
DZs)  and  at  approximately  60  percent  and  57  percent  for  Cases  B  and  C,  respectively  (in 
which  only  <2n£  SLM  has  DZ).  These  values  compare  quite  closely  to  those  of  33  percent  and 
52  percent  for  Cases  A  and  (B,  C).  respectively,  under  the  ODZ  condition  with  no  dc  block. 

The  6db  condition  cannot  be  designated  for  the  CMF  because  its  peaks  are  swamped  by  noise 
at  small  TDZ.  These  results  clearly  indicate  the  robust  operation  of  the  POF  and  BPOF 
while  operating  with  TDZ  noise. 


4.2  Correlation  Slgnal-To-Noise  Ratio 


The  definition  of  the  SNR  is 


The  numerator  is  the  ODZ  value  and  is  identical  in  value  to  that  in  Eq.  (1)  of  Reference  1. 
The  denominator,  however.  Includes  the  sum  of  all  four  k  terms  defined  in  Eqs.  (12)-(15). 

The  number  of  correlation  points  whose  amplitudes  (k-plj  are  less  than  one-half  of  (kp)jnax 

given  by  n-j-. 

The  SNR  simulation  results  for  the  three  cases  utilizing  the  three  filters  in  each  case  and 
a  1 -pixel  dc  block  in  the  filter  plane  are  shown  in  Figure  6.  Solid  curves  are  used  for  the 
POF  and  CMF.  while  dashed  curves  are  used  for  the  BPOF.  The  dotted  portions  of  these 
curves  indicate  that  a  secondary  peak  has  become  higher  in  amplitude  than  the  correlation 
peak.  10  The  vertical  bars  on  the  right  hand  side  of  Figure  6  indicate  the  range  in  variation 
in  the  SNR  simulations  using  the  three  filters  under  the  ODZ  conditions  when  no  dc  block 
was  used.i 

The  curves  in  Figure  6  show  that  SNR  falls  off  significantly  as  TDZ  increases.  This  is 
due  to  the  increasing  background  noise.  Such  fall-off  is  in  marked  contrast  to  the  relative 
flatness  of  the  SNR  plots  when  the  DZ  is  opaque,  as  indicated  by  the  vertical  bars  on  the 
light  side  of  the  figure.  For  a  fixed  value  of  TDZ  the  POF  is  the  best  performing  filter  and 
the  CMF  is  the  worst.  Also,  performance  is  more  adversely  affected  when  the  input  has  TDZ 
(as  in  Cases  A  and  B)  rather  than  the  filter  (as  in  Case  O.^*  The  substantial  energy  coming 
through  input  plane  TDZ  makes  a  large  noise  contribution  while  the  signal  is  decreased.  In 
fact,  for  input  TDZ  of  about  8  percent  the  noise  energy  from  the  (64x64)  pixel  TDZ  is  equal  to 
the  energy  in  the  (32x32)  pixel  input  signal. 
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TRANSMISSIVE  DEAD  ZONE,  Zi  (%) 


Figure  6.  Correlation  SNR  Versus  Transmissive  Dead  Zone  Area 
(in  Percent)  Using  Each  of  the  Three  Filters  for  Cases  A,  B  and  C. 
The  dotted  portions  of  each  curve  show  when  a  subsidiary  peak  is 
higher  in  amplitude  than  the  correlation  peak.  The  vertical  bars 
on  the  right  indicate  the  SNR  variations  for  each  filter  under 
opaque  dead  zone  conditions  and  no  dc  block. 


15 


L 


Table  2  demonstrates  how  the  dc  block  removes  a  considerable  portion  of  the  back¬ 
ground  noise  while  slightly  reducing  the  correlation  peak.  The  net  result  Is  a  significant 
Increase  In  the  SNR.  ^  *  In  the  central  portion  of  the  table  we  list  ratios  of  the  SNR  obtained 
with  a  1 -pixel  dc  block  in  the  system  to  the  SNR  obtained  with  no  block  under 
representative  conditions  of  TDZs.  Inclusion  of  this  dc  block  increases  the  SNRs  typically 
by  factors  of  1. 1  to  3.3.  However,  we  note  that  the  greatest  effect  occurs  for  Cases  A  and  B,  in 
which  SLMl  contains  TDZ,  rather  than  for  Case  C.  In  which  there  is  TDZ  only  In  SLM2. 

The  ratios  In  parentheses  pertain  to  those  conditions  when  a  secondary  peak  is  higher  than 
the  correlation  peak.  On  the  right  side  of  Table  2  we  list  the  ratios  of  correlation  peak 
Intensities  under  the  same  conditions  as  above,  namely,  with  and  without  a  dc  block  in  the 
system.  As  expected,  inclusion  of  a  dc  block  diminishes  Ip.  However,  these  reduced  peaks 

are  still  greater  than  the  true  correlation  peaks  given  by  k^.  We  would  also  note  that  the  dc 
block  reduction  of  the  TDZ  background  terms  of  Eqs.  (13)-(15)  is  equivalent  to  the 
Improvement  in  correlation  SNR  reported  by  Flannery,  et  al^^  and  the  reduction  in  clutter 
noise  observed  by  Fielding  and  Horner.*^ 


Table  2.  Comparison  of  SNR  and  Peak  Intensity  for  TDZ 
With  and  Without  a  1 -Pixel  dc  Block 


(SNR)rfeb/{SNR 

no  deb 

^lo^dcb/^lo^no  deb 

Case 

DZ  (%) 

POF 

BPOF 

CMF 

POF 

BPOF 

CMF 

A,B,C 

0 

1.2 

1.1 

1.7 

0.97 

0.99 

0.72 

A 

25 

2.3 

2.2 

(3.1) 

0.82 

1.00 

(0.23) 

B 

44 

2.1 

1.6 

3.3 

0.75 

0.94 

0.17 

C 

44 

1.3 

1.2 

(1.4) 

0.96 

1.02 

(0.70) 

Actually,  the  beam  focused  on  the  center  of  SLM2  has  a  half-power  beamwidth  of  about  3 
pixels.  Consequently,  a  dc  block  that  is  3  pixels  square  would  significantly  improve  the 
SNR,  even  though  It  would  lower  Ip  still  further. 

Similar  information  normalized  to  ODZ  values  (without  a  dc  block)  is  supplied  in  Table 
3.  Apart  from  the  values  at  0  percent  DZ  we  see  that  TDZ  in  any  SLM,  even  with  a  dc  block 
present,  diminishes  the  SNR  significantly  compared  to  that  for  the  ODZ  condition.  The 
right  side  of  the  table  shows  that  TDZ  enhances  the  peak  intensities  relative  to  the  ODZ 


Flannery,  D.L..  Loomis.  J.S.  and  Milkovlch,  M.E.  (1988)  Transform-ratio  ternary 
phase-amplitude  filter  formulation  for  improved  correlation  discrimination.  Appl.  Opt, 
27:4079-4083. 

12  Fielding,  K.H.  and  Homer  J.L.  (1990)  1-f  binary  joint  transform  correlator.  Opt  Eng., 
29(9):  1081-1087. 
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values  by  about  10  percent  for  the  POF,  by  10-30  percent  for  the  BPOF,  and  by  40-80  percent 
for  the  CMF. 


Table  3.  Comparison  of  SNR  and  Peak  Intensity  for  TDZ 
(With  1 -Pixel  dc  Block)  Relative  to  ODZ 


(SNR),,,b/(SNR)oDZ 

(ydcb/flo^ODZ 

Case 

DZ  (%) 

POF 

BPOF 

CMF 

POF 

BPOF 

CMF 

A.B.C 

0 

1.2 

1.1 

1.7 

0.97 

0.99 

0.72 

A 

25 

0.57 

0.47 

(0.42) 

1.12 

1.22 

(1.83) 

B 

44 

0.33 

0.31 

0.31 

1.15 

1.26 

1.72 

C 

44 

0.82 

0.75 

(0.56) 

1.06 

1.10 

(1.45) 

4.3  Energy  Throughput  Efficiency 


Table  4  shows  what  happens  to  the  optical  input  energy  at  different  locations  in  the 
system  as  TDZ  varies.  We  show  results  only  for  the  POF  with  a  1 -pixel  dc  block  in  the  filter 
plane  for  Cases  A.  B,  and  C.  Values  are  stated  in  terms  of  percentages  of  the  total  light 
exiting  from  SLMl.  The  values  in  parentheses  are  the  result  when  no  dc  block  is  present. 
The  values  in  the  last  four  columns  refer  to  the  following  four  quantities,  respectively:  light 
that  SLMl  diffracts  at  too  large  an  angle  so  that  it  misses  the  filter:  light  eliminated  by  the 
dc  block;  light  that  SLM2  diffracts  at  too  large  an  angle  so  that  it  misses  the  detector;  and 
light  recorded  by  the  detector.  This  latter  value  is  the  system's  energy  throughput  efficiency 
and  is  also  related  to  the  peak-to-correlation  energy  (PCE).^^  With  the  exception  of  the  light 
eliminated  by  the  dc  block,  the  other  three  partitions  of  energy  are  the  TDZ  counterparts  of 
curves  (2).  (4).  and  (5)  in  Figure  8  of  Reference  1.  The  values  across  any  row  always  add  to 
100  percent.  Note  that  Case  C  gives  the  highest  throughput  efficiencies. 


Kumar.  B.V.K.  and  Hassebrook.  L.  (1990)  Performance  measures  for  correlation  filters. 
AppL  Opt,  29:2997-3006. 
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Table  4.  Energy  Partitioning  (Percent)  Throughout  Correlator 
With  a  POP  and  a  1 -Pixel  dc  Block 


(Values  in  parentheses  are  for  no  dc  block) 


%  MISSING 
THE  FILTER 


.9 


(4.9) 


68.1 


64.4 


(64.4) 


50.0 


34.0 


23.0 


0.6 


.9 


(4.9) 


68.1 


64.4 


50.1 


(50.1) 


34.0 


23.0 


0.6 


4.9 


(4.9) 


4.9 


(4.9) 


4.9 


(4.9) 


.9 


(4.9) 


4.9 


(4.9) 


.9 


(4.9) 


4.9 


%  dc 
BLOCKED 

%  MISSING 
THE  DETECTOR 

13.2 

29.3 

(0) 

(29.2) 

16.5 

10.6 

19.5 

11.4 

(0) 

(12.6) 

28.8 

15.0 

38.2 

19.2 

44.6 

21.4 

33.2 

36.8 

%  RECORDED 
BY  DETECTOR 


29.3 


(29.2) 


8.4 


9.5 


13.9 


(19.9) 


18.9 


22.3 


29.2 


13.2 

29.3 

(0) 

(29.2) 

13.2 

45.4 

(0) 

(44.5) 

13.2 

45.8 

(0) 

(43.9) 

13.2 

45.9 

(0) 

(42.7) 

13.2 

41.7 

(0) 

(37.6) 

13.2 

33.5 

(0) 

(29.0) 

10.9 

5.6 

52.6 


(65.9) 


4.9 


4. 


(23.0) 


6.1 


8.5 


11.0 


29.4 


52.6 


(65.9) 


7.0 


6.5 


.2 


(30.0) 


8.8 


10.1 


12.9 


52.6 


(65.9) 


36.5 


(50.7) 


36.2 


(51.3) 


36.0 


(52.4) 


40.3 


(57.5) 


48.5 


(66.1) 


78.7 


5.  ANALYSIS  OF  CORRELATOR  PERFORMANCE 


The  qualitative  behavior  of  Ip  and  SNR  with  change  in  TDZ  can  be  demonstrated  from 
an  analysis  of  Eqs.  (12)-(15). 


5.1  Peak  Intensity 

The  total  correlation  output  amplitude  k-j-  is  given  by  the  sum  of  Eqs.  (12)-(15)  and  the 
corresponding  correlation  intensity  by  Ikjp.  Let  us  write  this  intensity  as  follows: 


0-^2)  \  \  \  +  Z,  (1-^2) 


(17) 


where  the  A's  represent  the  bracketed  terms  in  Eqs.  (12)-(15).  The  quantity  Ajj  is 
Independent  of  TDZ  amount  (as  well  as  of  the  filter  function)  because  it  depends  only  on  the 
q  and  r  functions  of  Eqs.  (8)  and  (10):  whereas,  the  other  three  A  quantities  have  a  slight 
dependence  on  Z  through  the  functions  v  and  m  of  Eqs.  (7)  and  (9).  Nonetheless,  to  a  first 
approximation,  we  can  consider  the  A’s  to  be  constants  with  respect  to  Z. 

We  want  to  show,  as  demonstrated  in  Figures  3-5,  that  peak  intensity  should  decrease 
as  TDZ  Increases.  As  an  example,  consider  Case  C.  Using  Eq.  (17)  let  us  take  the  ratio  of 
peak  intensity  Ip  for  eurbitrary  TDZ  Z2  to  Ip  at  Z2=0: 


(^2) 

lp(0) 


ri  -  A  ^ 

1  - 

s 

z. 

A 

2 

^  c  / 

(18) 


where  the  A's  are  computed  at  the  location  of  the  correlation  peak.  For  input  signals  of  the 
type  used  here  and  for  the  two  phase  filters.  is  usually  much  less  than  A^.  (at  least  over 
most  of  the  range  of  Z).  Consequently,  we  can  see  from  Eq.  (18)  that  as  Z2  increases,  the 
above  ratio  diminishes,  indicating  that  lp(Z2)  is  getting  progressively  smaller  than  lp(0). 

Similar  arguments  can  be  made  for  Cases  A  and  B.  thereby  confirming  the  behavior  of 
the  two  phase  filter  curves  in  Figures  3-5. 

For  the  CMF  the  relative  values  of  the  A  terms  are  sufficiently  different  from  those  of 
the  phase  filters  to  produce  a  minimum  in  the  Ip-Zj  curves  when  the  filter  SLM  contains 
TDZ.  In  Figures  3  and  5  we  see  that  these  minima  occur  in  the  vicinity  of  20-25  percent  DZ. 
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5.2  Signal'to-Noise  Ratio 


We  can  give  an  approximate  analysis  to  show  why  the  SNR  curves  behave  as  shown  in 
Figure  6.  Again,  we  take  Case  C  as  an  example.  Using  the  terms  of  Eq.  (17)  in  the  definition 
of  SNR  from  Eq.  (16)  and  dividing  numerator  and  denominator  by  (I-Z2),  we  get  for  this  case: 


SNR  (Case  C)  = 


(^)p 


nl/2 


(19) 


where  (A<,)p  is  the  value  of  A^.  at  the  location  of  the  correlation  peak  and  the  prime  indicates 
that  the  summation  is  to  include  only  those  points  in  the  correlation  pattern  whose 
amplitudes  are  less  than  one-half  of  the  correlation  peak.  Since  the  A^.  and  the  A^  have 
only  a  weak  dependence  on  the  amount  of  DZ,  as  Z2  increases,  the  denominator  increases 
while  the  numerator  remains  relatively  constant.  Hence,  regardless  of  the  filter  type,  we 
can  see  from  this  expression  that  SNR  diminishes  with  increase  in  TDZ. 

A  similar  type  of  argument  would  hold  for  the  SNRs  in  Cases  A  and  B. 

Figure  6  also  shows  that  for  each  filter  type  the  SNRs  for  the  three  cases  are  always  in 
the  order:  SNR(Case  C)>SNR(Case  B)>SNR(Case  A).  We  can  show  that  this  occurs  upon 
considering  the  ratios  of  various  SNRs.  For  excunple,  if  we  take  the  ratio  of  SNR(Case  C)  to 
SNR(Case  B),  with  both  of  them  expressed  in  the  form  suggested  by  Eq.  (19),  we  get,  after 
cancellation  of  the  common  terms 


SNR  (CASE  C) 
SNR  (CASE  B) 


,1/2 


L  J 


1/2 


(20) 


Let  us  consider  this  ratio  at  some  common  TDZ  so  that  Z2=Z2.  Then,  since  Ap^  is  greater 
than  Ag  for  the  types  of  images  used  here,  the  numerator  is  always  greater  than  the 
denominator.  Consequently,  SNR(Case  C)  is  always  greater  than  SNR(Case  B)  for  any 
TDZ>0. 

The  same  argument  can  be  used  to  show  that  SNR(Case  B)  is  always  greater  than 
SNR(Case  A)  for  any  TDZ>0. 
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5.3  Performance  With  a  dc  Block 


Equation  (A-20)  in  the  Appendix  shows  that  when  a  dc  block  is  included,  the  value  of 
each  A  term  in  Eq.  (17)  is  diminished  because  a  non-zero  term  has  been  subtracted  from  the 
original  value.  This  means  that  when  a  dc  block  is  included  all  of  the  A  terms  have  smaller 
values  than  they  would  have  without  the  block.  When  this  information  is  inserted  into  the 
appropriate  SNR  equation,  as  shown,  for  example,  in  Eq.  (19)  for  Case  C.  we  can  see  that  the 
denominator  will  decrease  faster  than  the  numerator.  Consequently,  the  ratio  will  have  a 
greater  value  than  it  would  have  had  without  the  dc  block. 


6.  DESIGNS  FOR  REDUCING  FALSE  PEAKS  AND  BACKGROUND 


The  presence  of  TDZ  in  the  input  plane  introduces  severe  noise  in  the  correlation 
pattern,  resulting  in  false  correlation  peeiks,  additional  background  patterns  and  a  SNR  that 
depends  on  the  correlator  design.  We  now  address  each  of  these  problems  in  turn. 


6. 1  Elimination  of  False  Correlation  Peaks 

When  the  reference  image  used  to  produce  the  filter  is  the  same  size  as  SLMl  and  this 
SLM  contains  TDZ,  then  false  peaks  eu'e  produced  in  the  correlation.  The  matching  of  their 
edges  contributes  significantly  to  the  overall  pattern,  as  well  as  adding  a  large  contribution 
to  the  correlation  peak.  As  an  example,  we  zmalyzed  a  Case  B  situation  in  which  SLMl  had 
44  percent  TDZ,  the  POF  in  the  filter  plane  was  the  same  size  as  SLMl  and  there  was  no  dc 
block.  Even  with  no  signal  in  the  input  plane  (s=0)  a  "correlation  peak"  was  produced  whose 
amplitude  was  just  as  large  as  that  from  a  correlation  between  the  POF  cmd  the  regular  input 
signal  having  no  TDZ.  Under  these  conditions,  the  total  correlation  amplitude  would  reduce 
to  fcp=kpfj.  We  found  that  such  false  contributions  to  the  correlation  peak  could  be 

eliminated  by  mismatching  the  sizes  of  the  reference  image  and  the  input  SLM.  We  finally 
made  the  reference  image  one-half  as  large  as  the  input  SLM  cmd  positioned  the  image  so 
that  its  edges  (32x32  pixels)  were  well  inside  of  the  SLMl  edges  (64x64  pixels). 


6.2  Background  Reduction 

Since  the  BPOF  correlation  consists  of  two  terms  it  is  instructive  to  write  its  k^.  term  as 
kc=lq,i+lq.2.  where  the  subscript  1  refers  to  the  correlation  portion  of  k^  and  the  subscript  2 
to  the  convolution  portion.  For  the  CMF  and  POF  correlators,  k^.2  is  zero. 
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If  the  centers  of  the  input  image  and  the  filter  reference  image  cire  located  on  the  optic 
axis  and  both  SLMs  have  TDZs,  then  all  of  the  terms  in  Eqs.  (12)-(15)  (including  kd 
and  kc2  for  a  BPOF)  will  also  have  their  centers  on-axis  in  the  correlation  plane. 
Consequently,  the  three  additional  noise  terms  of  Eqs.  (13)-(15).  plus  the  k^2  portion  for  a 
BPOF,  will  coherently  add  significant  contributions  onto  amd  around  the  correlation  peak. 
But,  if  the  filter  is  fabricated  from  a  filter  reference  image  whose  center  has  been  diagonally 
offset  from  the  optic  axis,  then  the  centers  of  k^.,  and  k^2  will  move  in  opposite  directions 
along  this  diagonal  while  the  centers  of  kg,  kj^  emd  kpj^  will  remain  fixed  on  the  optic  axis. 
Moreover,  if  the  input  signal  is  translated,  then  the  combination  of  k^i,  k(.2  and  kg  will 
translate  as  one  fixed  pattern.  This  process  provides  us  with  a  means  of  separating  the  peak 
of  kci  from  the  background  noise  of  k^2  mid  kg.  The  same  separation  phenomenon  was  also 

observed  in  amplitude-encoded  BPOFs.^'^-*^ 

In  our  case,  although  the  entire  field  in  each  plane  extends  to  128x128  pixels,  the  actual 
dimensions  of  the  signal,  the  filter  and  the  detector  cover  64x64  pixels.  The  patterns  of  kg, 
kR,  and  kpR  cover  32x32,  64x64  and  somewhat  greater  than  64x64  pixels,  respectively.  We 
shifted  the  filter  reference  image  into  a  corner  quadrant  in  the  input  plane  by  a  distance  of 
32  pixels  along  the  x  and  y  directions.  For  the  BPOF  this  reduced  the  backgrounds  from 
Eqs.  (12)  and  (13)  since  kj.i  is  shifted  to  a  low  intensity  comer  of  both  the  kg  image  and  the 
kc2  convolution. 

Additional  noise  reduction  is  provided  by  the  dc  block,  which  reduces  the  low  frequency 
energy  of  the  overlapping  kp  and  kpR  (as  well  as  the  kgand  kp2)  contributions.  Consequently, 
both  the  shifting  process  and  the  dc  block  allow  the  peak  to  become  relatively  insensitive  to 
the  noise  terms. 

Downie  and  Reid*^  have  reported  that  for  their  BPOF  correlator,  the  correlation  peak 
remained  relatively  independent  of  changes  in  threshold  line  angle.  We  studied  the  effects 
of  changing  the  BPOF  threshold  line  angle  in  Case  A  with  Zi=Z2=44  percent  and  found  a  +17 
percent  maximum  variation  in  SNR,  as  well  as  a  +25  percent  variation  in  Ip  about  the  tme 
correlation  value.  In  Figure  3  this  would  mean  that  the  value  on  the  dashed  curve  at  44 
percent  TDZ  could  vary  between  ±25  percent  around  the  value  on  the  solid  curve.  The  Ip 
fluctuation  came  from  coherent  interference  between  the  correlation  peak  and  the 
background  terms  as  their  relative  phases  varied  with  threshold  angle.  The  SNR  fluctuation 
came  from  the  phase  variations  of  kp2  and  kpR  as  the  change  in  threshold  angle  shifts  the 
Intensity  distribution  in  the  correlation  plane  relative  to  our  64x64  pixel  detector.  The 
substantial  improvement  in  SNR  (together  with  a  decrease  in  Ip  for  a  centered  filter 
reference  image)  observed  by  Downie  and  Reid’^  results  from  the  appearance  and 
disappearance  of  a  dc  block  in  the  amplitude-encoded  BPOF  as  the  threshold  angle 


14  Flavin.  M.A.  and  Homer.  J.L.  (1988)  Amplitude  encoded  binary  phase-only  filters, 

Proc.  S.P.I.E.,  938:261-265,  Digital  and  Optical  Shape  Representation  and  Pattern  Recognitioru 
R  Juday  (ed.). 

15  Downie,  J.D.  and  Reid,  M.B.  (1990)  Mapping  considerations  for  optimal  binary  correla¬ 
tion  filters,  Appl.  Opt,  29:5235-5241. 


22 


changes.^®  Flanneiy,  et  al, have  also  reported  similar  variations  for  the  BPOF  correlation 
of  objects  in  clutter. 

A  3D  correlation  pattern  for  our  system  using  a  reference  image  mismatched  to  the  SLM 
size,  as  well  as  off-axis  shifts  for  the  reference  image  and  input  signal,  as  described  above,  is 
shown  in  Figure  7.  Here,  the  system  has  a  BPOF,  a  1 -pixel  dc  block  and  no  TDZ  in  SLM2 
and  a  TDZ  of  44  percent  in  SLM  1  (Case  B).  The  ridges  surrounding  the  peak  are  due  to  the 
edges  of  the  dc  filtered  kp^  term  plus  the  correlation  of  the  64x64  pixel  SLMl  edges  with  the 

32x32  pixel  reference  image  edges. 


Figure  7.  A  Representative  3D  Pattern  for  a  Correlator  Having 
a  BPOF,  a  1 -pixel  dc  Block  With  No  TDZ  in  SLM2  and  a  TDZ  of 
44  Percent  in  SLM  1  (Case  B).  The  reference  image  is  one-half  as 
large  as,  and  is  centered  within,  the  input  SLM.  The  SNR  is  6.4. 


6.3  Design-Dependent  SNR 


The  quantitative  value  that  one  obtains  for  the  SNR  depends  on  how  one  chooses  the 
spatial  region  over  which  the  rms  noise  in  the  denominator  is  calculated.  We  have  followed 
a  standard  convention  of  choosing  a  64  square  pixel  input,  filter  and  detector.  In  our  design 
we  chose  the  detector  position  by  centering  it  on  the  correlation  peak  when  the  input  image 
was  centered  in  SLMl.  With  this  design  the  peak  will  stay  on  the  detector  as  long  as  the 


1®  Khouiy,  J.,  Gianino,  P.D.,  Kane,  J.  and  Woods,  C.L.  (1994)  Edge  enhancement  techniques 
for  improving  the  performance  of  binary  phase-only  filter  pattern  recognition  devices.  Opt. 
Eng.,  33(31:856-864. 

Flannery,  D.L.,  Loomis,  J.S.  and  Milkovich,  M.E.  (1988)  Design  elements  of  binary 
phase-only  correlation  filters,  AppL  Opt.,  27:4231-4235. 
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center  of  the  input  Image  is  within  SLM 1 .  We  use  this  detector  area  as  the  noise  area  for  the 
rms  calculation. 

Our  detector  array  contains  most  of  the  (0.  0)  order  of  the  correlation  pattern,  essentially 
all  of  the  kg  image,  (approximately  1/4  of  k^.2  for  the  BPOF)  and  roughly  1/2  of  the  kp  and 
kpR  terms.  It  also  contains  parts  of  the  neighboring  higher  orders  of  the  correlation 
pattern.  The  center  portion  of  the  area  shown  in  Figure  7  corresponds  to  our  detector,  while 
the  entire  area  is  twice  that  of  our  detector.  Furthermore,  the  reference  image  shift  and 
input  image  position  will  determine  what  portions  of  the  four  k  terms  will  fall  within  the 
detector.  Clearly  then,  one  can  see  that  detector  size  affects  the  calculated  value  of  SNR. 


7.  MIXED  CASES  OF  TDZ  AND  ODZ 


It  is  also  possible  to  model  mixed  cases  analytically,  such  as  when  one  SLM  contains 
transmissive  dead  zone  (TDZ)  only  and  the  other  opaque  dead  zone  (ODZ)  only  or,  when  any 
one  SLM  contains  both  TDZ  and  ODZ.  Consider,  for  example,  the  first  situation  in  which 
the  input  SLM  has  ODZ  of  Zj  and  the  filter  SLM  has  TDZ  of  Z2.  We  call  this  the  (ODZ.TDZ) 

case.  Here,  tglxj.  yj)  =  t^  (X3,  yg)  =  0.  So,  from  Eqs.  (3)-(6)  we  see  that  only  the  k^  and  kg 
terms  exist.  Upon  applying  this  information  to  Eq.  (17)  we  can  simplify  that  equation 
(referring  to  Table  1)  to: 


|k,j,|  VdZ.  TDZ)  =  (1  -  zf  |(  1  -Z^)  +  Z^AJ^  =  (1  -  zf  •  [|kj.|  Vase  C)]  . 


(21) 


Hence,  the  peak  intensities  for  this  situation  can  be  obtained  by  multiplying  the  data  in 
Figure  5  by  the  factor  (1-Zi)2.  Moreover.  SNR  for  this  situation  is  identically  equal  to  the 
SNR(Case  C)  from  Figure  6  because  multiplication  of  k]  (Case  C)  by  the  constant  value  (l-Zj) 
does  not  affect  the  SNR. 

In  the  reverse  situation  in  which  the  input  SLM  has  TDZ  of  Zj  and  the  filter  has  ODZ  of 
Z2.  which  we  call  the  (TDZ.ODZ)  case,  then  Tp(^2-  ^2)  =  ya)  =  0-  From  Eqs.  (3)-(6)  we  see 
that  only  the  k^.  and  kpp  terms  exist.  Applying  this  information  to  Eq.  (17)  we  can  simplify 

that  equation  to 


|k.j,|  VdZ.  ODZ)  =  (1  -  Z^)^  |(  1  -  Zj)  A^  +  ZjAj^J^  =  (1  -  Z^f  [k^l  Vase  B)]  . 


(22) 
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Consequently,  the  peak  intensities  for  this  situation  can  be  calculated  by  multiplying  the 
data  in  Figure  4  (Case  B)  by  the  factor  (1-22)^  and  the  SNR  is  identically  equal  to  SNR(Case 

B). 

For  the  (ODZ.TDZ)  case  no  dc  block  is  necessary  and  no  false  correlation  peaks  due  to 
edge  correlation  (see  Section  4)  occur;  whereas,  for  the  (TDZ,ODZ)  case  a  dc  block  is 
necessary  and  contributions  from  edge  correlation  can  occur.  So,  if  one  has  a  choice 
between  these  two  combinations,  it  is  generally  better  to  place  the  ODZ  SLM  in  the  input 
plane. 

The  most  general  situation  occurs  when  the  DZ  of  any  one  SLM  consists  of  both  TDZ  and 
ODZ,  that  is,  (DZ)j=(TDZ)i+(ODZ)i.  In  fact,  all  of  the  preceding  analyses  can  be  regarded  as 
special  cases  of  this  situation,  with  either  (ODZ)j  or  (TDZ);  equal  to  zero.  This 
generalization  can  be  very  easily  accounted  for  by  considering  Eqs.  (12)-(15).  There,  the  only 
terms  that  Involve  DZ  directly  are  the  Z;  and  (l-Zj),  which,  of  course,  refer  to  (TDZ);  and 
(AZ)j,  respectively.  So  we  see  that  the  only  portions  of  DZ  that  affect  these  equations  are  the 
amounts  of  TDZ  and  AZ;  the  opaque  amount  of  DZ  plays  no  role  here  at  all.  Hence,  if  in  all 
of  the  previous  equations  we  interpret  all  Z;  to  refer  to  TDZ  only  and  all  (1-Z;)  to  refer  to  AZ 
only,  then  all  of  the  previous  equations  will  be  applicable  to  this  general  case.  Note  that  the 
kj,  still  represents  the  "true"  correlation  because  it  depends  on  the  amounts  of  AZs  only  and 
is  Independent  of  the  amounts  of  TDZ  and/ or  ODZ  that  may  be  present. 


8.  CONCLUSIONS 


Analytical  results  show  that  there  are  four  basic  terms  making  up  the  total  correlation: 
one  gives  the  true  correlation,  which  is  the  same  as  that  for  the  corresponding  ODZ  mode, 
and  the  other  three  form  the  background  noise. 

When  the  results  from  the  mathematical  analysis  are  combined  with  those  from  the 
computer  simulation,  we  show  that  for  a  fixed  value  of  TDZ: 


1.  Correlation  performemce  in  the  form  of  SNR,  peak  intensity  (Ip)  and  throughput 
efficiency  is  always  the  highest  when  a  POF  is  used  in  the  system,  is  slightly  poorer 
for  a  BPOF,  and  is  very  poor  for  a  CMF. 

2.  Insertion  of  a  1 -pixel  dc  block  at  the  center  of  the  filter  plane  SLM  enhances  the  SNR 
by  removing  the  background  in  the  correlation  plane.  The  enhancement  factors  are 
from  2  to  3  when  the  TDZ  is  in  the  input  SLM  and  about  1.5  when  there  is  no  DZ  at 
all  or  when  the  TDZ  is  in  the  filter  plane  SLM  only. 
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3.  For  either  of  the  phase  filters,  Ip(Case  B)>lp  (Case  C)  >  Ip(Case  A):  whereas,  for  all 
three  filters,  SNR(Case  C)  >  SNR(Case  B)  >  SNR(Case  A). 

4.  For  either  of  the  phase  filters  the  Ip  for  the  TDZ  mode  are  always  greater  than  those 
for  the  ODZ  mode  even  when  a  1 -pixel  dc  block  is  used.  This  is  not  so  for  the  CMF. 

5.  For  an  ODZ  input  and  a  TDZ  filter  the  correlation  is  reduced  from  that  of  Case  C  by  a 
factor  of  (1-Zi)2,  but  the  SNR  is  the  same  as  in  Case  C. 


6.  For  a  TDZ  input  and  an  ODZ  filter  the  correlation  is  reduced  from  that  of  Case  B  by  a 
factor  of  (1-Z2)2.  but  the  SNR  is  the  same  as  in  Case  B. 

7.  For  an  SLM  containing  both  ODZ  and  TDZ  the  mathematical  analysis  derived  here 
also  holds,  except  that  the  term  Z,  pertains  only  to  the  amount  of  TDZ  in  each  pixel 
and  the  term  (l-Zj)  pertains  only  to  the  amount  of  AZ. 

8.  The  efficiency  of  the  signal  light  getting  through  to  the  detector  is  much  greater  for 
Case  C  than  it  is  for  Cases  A  or  B. 


As  TDZ  Increases,  the  three  salient  results  coming  from  the  simulation  are: 


1.  SNRs  decrease  by  almost  an  order  of  magnitude.  This  is  completely  different  from 
the  ODZ  mode  in  which  the  SNR  was  not  affected  at  all  by  changes  in  ODZ  area. 

2.  The  correlation  peak  intensities  Ip,  although  increased  by  TDZ  noise,  drop  off  by 
factors  of  10  to  100,  just  as  they  did  in  the  ODZ  case.  And  these  intensities  vary 
slightly  with  changes  in  filter  input  position,  as  they  did  in  the  ODZ  case. 

3.  For  sufficiently  large  TDZ,  secondary  peaks  in  the  correlation  pattern  become  larger 
than  Ip  and  keep  increasing  in  amplitude  up  to  100  percent  DZ.  This  is  due  to  the 
large  amount  of  light  passing  through  the  DZs.  For  Case  A  this  occurs  after 
approximately  60  percent  DZ  for  both  of  the  phase  filters  and  after  approximately 
20  percent  for  the  CMF. 
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One  can  always  use  the  following  as  an  approximate  rule-of-thumb:  Ip  drops  off  by  6db 
(a  factor  of  4)  below  its  zero  DZ  values  when  only  one  SLM  has  an  active  zone  in  the  vicinity 
of  40  percent  or  when  both  SLMs  have  an  active  zone  of  61  percent  each. 


9.  RECOMMENDATIONS  TO  USERS 

1.  One  should  always  try  to  place  the  SLM  with  the  smaller  amount  of  TDZ  in  the 
input  plane. 

2.  When  the  input  SLM  has  TDZ  it  is  advantageous  to  insert  a  dc  block  in  the  filter 
plane. 

3.  For  a  square  reference  image  it  is  absolutely  necessary  to  mismatch  the  reference 
image  and  SLM  sizes  to  minimize  the  false  correlation  enhancement  from 

the  matching  of  the  input  and  TDZ  edges. 

4.  To  allow  the  correlation  peak  to  become  relatively  insensitive  to  the  TDZ  (and 
possibly  convolution)  noise  terms,  one  should  fabricate  the  filter  from  a  filter 
reference  image  whose  center  has  been  offset  from  the  optic  cixis. 
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Appendix  A 

Derivation  of  Signals  w',  g,  tp  and  tjc 


Al.  NO  dc  BLOCK  IN  FILTER  PLANE 


The  signal  amplitude  w(xj,  yj)  encoded  on  SLMl  due  to  the  active  areas  in  its  pixels  is 
given  by:® 


w(Xi.yi)  = 


[{ 


I  I  8  (Xj  -  n^cj  5  (y^  -  n  c  )  s  (x^.  y^)  ^  O  (b^v^.  b  v^) 


n^=-oo  ny= 


} 


r  r 

X  y 


(A-1) 


where  s  is  the  band-limited  continuous  signal  that  would  produce  digital  samples.  To  model 
the  finite  size  of  the  active  area  of  the  SLM  pixels  the  sampled  s  is  convolved  with  the  2D 
rect  function  v^^Vy  of  Eq.  (7).  Then,  this  entire  function  is  multiplied  by  the  2D  rect  function 
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Txty  of  Eq.  (10),  which  accounts  for  the  physical  limits  (Lj^.  Ly)  of  SLMl.  This  result  is  the 
same  as  Eq.  (2)  in  Reference  1  except  that  the  full  2D  form  is  given  and  the  last  two  terms 
indicate  the  SLM's  physical  limitation. 

Before  continuing  with  this  analysis,  we  adopt  the  following  notation  for  certain 
functions  that  occur  often  in  the  filter  plane: 


=  sine  (bx  ^2)  =  ^ 

(A-2) 

9^  =  rect  (Cx  ^2)  =  ^ 

(A-3) 

=  rect  (hx  Lx  ^a/Cx)  =  (Cx/^x  ^  <nixl 

(A-4) 

R^  s  Lx  sine  (Lx  ^2)  =  {^xl 

(A-5) 

where  all  of  the  quantities  on  the  right  side  have  already  been  defined  and  the  Fourier 
transform  (FT)  equivalents  are  indicated.  There  will  be  another  identical  set.  namely,  V^, 
9^.  M^,  and  R^.  in  which  the  subscript  x  and  variable  ^2  are  replaced  by  y  and  ri2- 
respectively.  The  additional  terms  h^L^  and  hyLy  in  the  M  functions  of  Eq.  (A-4)  come  from 
the  factors  that  scale  space  coordinates  to  spatial  frequency  coordinates  (namely. 
hxLx^2  =  ^2  and  hyLyq2  =  y2)- 

The  Fourier  transform  of  w  in  Eq.  (A-1)  is  represented  by  W(^2'  ^2^-  model  SLMl  so 
that  all  orders  of  diffraction  (indexed  by  the  integers  n^.  Oy)  are  discarded  except  the  (0.  0) 
order  because  they  fall  outside  of  SLM2.  The  remaining  transform  is  given  as 


v.v 

^  n 


)OR,R 


Q.9  • 


(A-6) 


where  S  is  the  Fourier  transform  of  s  and  the  Q  terms  account  for  the  finite  extent  of  SLM2, 
whose  endpoints  occur  at  ^2  =  ±l/2Cj.  and  ^2  =  ±l/2Cy.  (When  the  term  is  assumed  to  be 
flat,  Eq.  (A-6)  gives  the  2D  form  of  Eq.  (4)  in  Reference  1).  Taking  an  inverse  Fourier 
transform  gives  the  amplitude  function: 


ya)  =  ( 1  -  2i)  [{ s  (X3.  ya)  O  Vy }  • 


r  r 
X  y. 


(A-7) 
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The  r  terms  show  that  w'  is  limited  by  the  physical  dimensions  L^,  Ly  of  the  detector. 

As  can  be  seen  from  Figure  1,  the  DZ  in  each  pixel  of  SLMl  can  be  represented  by  the 
continuous  function 


.  rect 

[''xJ 

c 

^  yJ 

-  b  V  .  b  V  . 
XX  y  y 


(A-8) 


Following  the  same  reasoning  as  in  Eq.  (A-1)  we  obtain  t^,  the  signal  emanating  from  SLMl 
that  comes  solely  from  the  TDZ,  as 


1 

M  8 

i  s(^-Vx)Myi-Vy) 

j-OD,(x,.y,) 

.  r  r  . 

X  y 

(A-9) 

L  L"x“-“ 

ny=-oo 

I  J 

Its  Fourier  transform,  Tg.  can  be  shortened  and  modified  to 


<  fe  "2)  =  V,  ■  SA  .  IA-1 

in  which  all  diffracted  orders  higher  than  (0,  0)  are  discarded  because  they  miss  the  filter, 
emd  the  physical  dimensions  of  SLM2  cU'e  accounted  for  in  the  Q  terms.  Taking  an  inverse 
Fourier  transform  yields  t^,  the  eimplitude  in  the  output  plane: 


lA-U) 


The  amplitude  of  the  filter  transform  function  G  1^2’  ^2^  encoded  on  SLM2  due  to  the 
active  areas  in  the  pixels  is  given  by: 
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G(^-^2)  = 


ny=-~ 


O  M,M 
4  n 


•  Q.Q 


(A- 12) 


where  F*  is  the  band-limited  continuous  filter  transform  function  for  any  one  of  the  three 
desired  filters  and  was  originally  derived  from  the  signal  w.  The  5  functions  sample  F*  at 
the  center  of  each  pixel  at  spatial  frequency  intervals  of  l/L^,  1/Ly.  These  intervals 
correspond  to  physical  widths  of  h,.,  hy,  respectively.  To  model  the  finite  size  of  the  active 
zone  of  the  SLM  pixels  the  sampled  F*  is  convolved  with  the  2D  rect  function  of  Eq. 

(A-4).  Then,  this  entire  function  is  multiplied  by  the  2D  rect  function  Q^Q^i  of  Eq.  (A-3) 
which  accounts  for  the  physical  limits  of  SLM2.  Equation  (A- 12)  is  the  2D  counterpart  of 
Eq.  (5)  of  Reference  1.  Taking  the  inverse  Fourier  transform  of  Eq.  (A- 12)  and  keeping  all 
diffraction  orders  for  completeness,  there  results 


(1-Z2) 


i 

.v=-» 


n  L  -y, 

XX  •'I 


+  n  L  ') 
y  yy 


m  m 

X  y 


} 


O  q  q 


(A- 13) 


where  f*(-x,-y)  is  the  inverse  Fourier  transform  of  F*(^,  q). 
The  DZ  in  each  pixel  of  SLM2  can  be  represented  by 


1^2  (^2- ^2)  =  ^^^^ 


rect 

.  rect 

-  rect  1-^  i .  rect 

h 

h 

1  e  j 

e 

\  xJ 

^  yJ 

V  X/ 

^  yJ 

(A- 14) 


or.  in  spatial  frequency  coordinates. 


^^2  (^’  ^2)  =  (^x^  •  ■ 


(A- 15) 
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The  signal  transmitted  by  the  TDZs  in  SLM2  is  given  by 


Tp  (^.  112)  " 


Q.Q 


(A- 16) 


Taking  the  inverse  Fourier  trcmsform,  evaluating  the  resulting  convolution  and  separating 
the  higher  diffraction  orders  from  the  zero  order,  yields 


yn  e 

sine 

X  X 

.  sine 

y  y 

h  . 

i  h 

V  X  J 

1  y  y 

sine 


n  L 

X  X 


sme 


y.  -  n  L 
•^3  y  y 


(A- 17) 


where  the  primes  on  the  summations  in  the  last  term  signify  that  n,.  and  n^  can  go  from  -«> 
to  +  00,  but  cannot  be  zero  simultaneously. 


A2.  WITH  A  dc  BLOCK  IN  FILTER  PLANE 


If  we  insert  a  rectangular  opaque  dc  block  of  p,.  by  Py  pixels  at  the  center  of  SLM2  (with 
dimensions  p^h^  by  pyhy),  its  transmittance  would  be 


rect 


(  y-? 


p  ^  , 
y  yy 


(A- 18) 


V 

This  would  scale  in  spatial  frequency  coordinates  to 


.  rect 

1  Px  J 

[Py  J 

(A- 19) 
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which  would  transform  to 


P.P 

‘<ic(’s-y3)=s(vy3)-in5®'"‘= 

X  y 


('Px’S 


sine 


Pyy3 


V  y  J 


(A-20) 


The  signals  G  and  Tp.  given  by  Eqs.  (A- 12  and  A- 16).  would  have  to  be  multiplied  by 
Tdc(^2’  ^2)-  This  would  result  in  rectangular  sections  being  blocked  out  of  their  centers.  The 
G  function,  for  example,  when  modulated  by  this  dc  block,  would  transform  as: 


G  ^2)  •  '^dc  (^-  ^2)  g  (^3-  ya)  ®  ‘dc  (^3-  ^3)  • 


(A-21) 


A  similar  expression  would  hold  for  Tp. 

This  dc  block  eliminates  all  the  lower  frequencies  of  light  in  the  ranges 


ly  -  -  Py/^^y  • 
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